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We present improved measurements of the branching fraction and CP asymmetry for the process 
p^n''\ The data sample corresponding to 211 fb~^ comprises 232 million T(4S) — + BB decays 
collected with the BABAR detector at the PEP-II asymmetric B Factory at SLAC. The yield and 
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CP asymmetry are measured using an extended maximum likelihood fitting method. The branching 
fraction and CP asymmetry are found to be B{B^ p^-k'^) = [10.2 ±1.4 (stat) ±0.9 (syst)] x 10"*^ 
and Acp{B^ ^ P*7r°) = -0.01 ±0.13 (stat) ±0.02 (syst). 

PACS numbers: 11.30.Er, 13.25.Hw 



Branching fraction and CP asymmetry measurements 
of charmless B meson decays provide valuable con- 
straints for the determination of the unitarity triangle 
constructed from elements of the Cabibbo-Kobayashi- 
Maskawa quark-mixing matrix P, In particular, the 
angle a = arg [—V^^V^^lV^^^*^^] of the Unitarity Triangle 
can be extracted from decays of the B meson to p^ir^ fi- 
nal states Q. However, the extraction is complicated by 
the interference of decay amplitudes with differing weak 
and strong phases. One strategy to overcome this prob- 
lem is to perform an SU(2) analysis that uses all pi: final 
states Assuming isospin symmetry, the angle a can 
be determined free of hadronic uncertainties from a pen- 
tagon relation formed in the complex plane by the five 
B pn decay amplitudes B — > p'^ir' - B° ^ 
B° p°7rO, B+ p+TT° and B+ p°7r+. These am- 
plitudes can be determined from measurements of the 
corresponding decay rates and CP asymmetries. While 
all these modes have been measured [H, Q, the current 
experimental uncertainties need to be reduced substan- 
tially for a determination of a. Here we present an update 
to previous measurements of the B^ — > p^Tr*^ branching 
fraction and CP asymmetry 

_ N{B- p-n") - N{B+ p+TT°) 
~ N{B- -> p-TT^) + N{B+ p+ttO) ■ 

The main additions compared to our previous anal- 
ysis [1| are a larger dataset, a study of possible back- 
grounds from higher p resonances and the use of the p 
mass in the maximum likelihood fit. 

The data were collected with the BaBAR detector Q 
at the PEP-H asymmetric-energy e'^e^ storage ring 
at SLAG. Charged-particle trajectories are measured 
by a five-layer double-sided silicon vertex tracker and 
a 40-layer drift chamber located within a 1.5-T mag- 
netic field. Charged hadrons are identified by combin- 
ing energy- loss information from tracking (dE/dx) with 
the measurements from a ring-imaging Cherenkov detec- 
tor. Photons are detected by a CsI(Tl) crystal elec- 
tromagnetic calorimeter with an energy resolution of 
aE/E 0.023(£;/ GeV)-^/'* © 0.014. The magnetic flux 
return is instrumented for muon and identification. 
The data sample includes 232 ± 3 million BB pairs col- 
lected at the T{AS) resonance, corresponding to an in- 
tegrated luminosity of 211 fb^^. In addition, 22 fb^^ of 
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data collected 40 MeV below the T(4S') resonance mass 
are used for background studies. We perform full detec- 
tor Monte Carlo (MC) simulations equivalent to 460 fb~^ 
of generic BB decays and 140 fb~^ of continuum quark- 
antiquark events (e"*"e~ — > gg, q = u,d, s, c). In addition, 
we simulate over 50 exclusive charmless B meson decay 
modes, including 1.4 million signal B^ — > p^Tr" decays. 

B meson candidates are reconstructed from one 
charged track and two neutral pions. The charged track 
used to form the B^ —> p^ir^ candidate is required to 
have at least 12 hits in the drift chamber, to have a 
transverse momentum greater than 0.1 GeV/c, and to be 
consistent with originating from the beam-spot. It must 
have ionization-energy loss and Cherenkov angle signa- 
tures consistent with those expected for a pion. We re- 
move charged tracks that pass electron selection crite- 
ria based on dE/dx and calorimeter information. Neu- 
tral pion candidates are formed from two photon can- 
didates, each with a minimum energy of 0.03 GeV and 
which are required to exhibit a lateral profile of energy 
deposition in the electromagnetic calorimeter consistent 
with an electromagnetic shower jT,]. The angular accep- 
tance of photon candidates is restricted to exclude parts 
of the calorimeter where showers are not fully contained. 
We require the photon clusters forming the tt" to be sep- 
arated in space, with a 7r° energy of at least 0.2 GeV and 
an invariant mass between 0.10 and 0.16GeV/c^. 

Two kinematic variables, AE = Eg — y/s/2 and the 
beam-energy substituted mass of the B meson toes = 

[s/2 + po • pb)^/Eq — p^, are used for the final se- 
lection of events. Here E'^ is the energy of the B me- 
son candidate in the center-of-mass frame, Eq and ^/s 
are the total energies of the e"'"e~ system in the lab- 
oratory and center-of-mass frames, respectively;and po 
and pb are the three-momenta of the e~^e~ system and 
the B meson candidate in the laboratory frame, re- 
spectively. For correctly reconstructed p^Tt" candidates 
AE peaks at zero, while for final states with a charged 
kaon, such as B^ K*^n^, AE is shifted by approx- 
imately 80 MeV on average. Events are selected with 
5.20 < TOES < 5.29GeV/c2 and |A£;| < 0.20 GeV. The 
AE limits remove background from two- and four-body 
B meson decays with a small loss in signal efficiency. 

Continuum events are the dominant background. To 
suppress this background, we select only those events 
where the angle ^fph in the center-of-mass frame between 
the sphericity axis 8] of the B meson candidate's decays 
products and the sphericity axis of the rest of the event 
satisfies jcosSgpjJ < 0.9. In addition, we construct a 
non-linear discriminant, implemented as an artificial neu- 
ral network (Ann) that uses three input parameters: the 
zeroth- and second-order Legendre event shape polyno- 
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mials Lq and L2 calculated from the momenta and polar 
angles, with respect to the B meson thrust axis, of all 
charged particle and photon candidates not associated 
with the B meson candidate, and the output of a mul- 
tivariate, non-linear B meson candidate flavor tagging 
algorithm The output Ann of the artificial neural 
network peaks at 0.5 for continuum- like events and at 
1.0 for B meson decays. We require Ann > 0.63 which 
reduces the continuum background by half for a 5% loss 
in signal MC efflciency. To further improve the signal- 
to-background ratio we restrict the invariant mass of the 
p candidate to 0.55 < 771^^-^ < 0.95 GeV/c^. 

The average B meson candidate multiplicity per event 
is 1.8 as neutral and charged pion combinatorics can lead 
to more than one B meson candidate. We choose the 
best candidate based on a formed from the measured 
masses of the two 7r° candidates within the event com- 
pared to the known 7r° mass • In the case of multiple 
charged pion candidates the choice is random so as not 
to bias the fit distributions. This random selection has a 
negligible impact on the systematic uncertainty. The to- 
tal B^ — > p^tt'^ selection efficiency is 15.4 ±0.1%. In sig- 
nal MC studies, the candidate is correctly reconstructed 
54.9% of the time. The remaining candidates come from 
self-cross-feed (SCF, 37.5%) and mistag events (7.6%). 
SCF events stem primarily from swapping the low en- 
ergy tt" from the resonance with another from the rest 
of the event. Signal events reconstructed with the wrong 
charge are classified as mistag events. Both SCF and 
mistag events emulate signal events, however the resolu- 
tion in rTiES and AE tends to be worse. 

We use MC events to study the backgrounds from other 
_B-meson decays. The dominant contribution comes from 
6 — > c transitions; the next most important is from 
charmless _B-meson decays. Seventeen individual charm- 
less modes show a significant contribution once the event 
selection has been applied. These modes are added into 
the fit (described below) fixed at the yield and asymme- 
try determined by the simulation, based on their mea- 
sured values [13]. The largest contributions come from 
B° ^ and B° p^Tr^. For B° r]'TT° and 

B^ K*'^p^ we use half the measured upper limit p^ . 
We estimate the B° o^tt'^ branching fraction from 
that of B^ a^Tr~ [ll| using isospin relations. If no 
charge asymmetry measurement is available, we assume 
zero asymmetry. 

Although all other states that decay like the p to tttt^ 
- the p(1450) and the /5(1700), subsequently referred to 
collectively as p* - lie outside our p(770) mass cut, a 
contribution to our signal cannot be ruled out a priori. 
To account for the possible presence of these modes, an 
unbinncd maximum likelihood fit to the B^ — > p*^n^ 
yield is performed in a sideband of the m^Tr invariant 
mass. This fit uses the same algorithm as described be- 
low but with only the three input variables mES, 
and Ann- The mass window is chosen to be as far as 
possible from the p(770) mass, centered near the pole 
of the p(1700) at 1.5 < m^^ < 2.0 GeV/c^. The fitted 



yield for the B^ p*^iT^ decay is then extrapolated 
into the p(770) region, 0.55 < m^^ < 0.95GeV/c^, using 
a non-relativistic Breit-Wigner line-shape. Although the 
choice of mass range is motivated by the p(1700), any 
yield seen is attributed entirely to the p(1450), which is 
the closer of the two resonances to the signal. From the 
B^ p*(1450)7r" MC, the ratio of the number of can- 
didates in the sideband to candidates in the signal mass 
region is approximately 12.6:1. The fit in the sideband 
yields 101 ± 32 events, resulting in an estimate of the p* 
background of 8 events. We investigate possible inter- 
ference effects by using an analytical model for the line- 
shapes of the p(770) and the p* . We compare the use of 
relativistic and non-relativistic Breit Wigner lineshapes 
and vary the widths of the lineshapes by their uncertain- 
ties [l3l • We also scan the relative phase between the two 
resonances from — tt to n. We assign a conservative sys- 
tematic uncertainty of 100% for the p* background based 
on the largest change in the number of events in the range 
0.55 < TO^^ < 0.95 GeV/c^ from these tests. The p* then 
enters into the nominal fit with PDFs constructed from 
B^ p±(1450)7rO MC simulation. 

The non-resonant B^ — s- tt^tt^tt^^ branching frac- 
tion has, to date, not been measured. To estimate 
the size of its contribution we select a region of the 
Dalitz plot — defined by the triangle (rn^±^o, ™^±^o) = 

(6, 6), (6, 15), (11, 11) GeV7c4 — that is far from the' sig- 
nal as well as the p(1450) and higher resonances and 
which has low levels of continuum background. The un- 
binned maximum likelihood fit with only three input vari- 
ables (toes I AE, and Ann) is applied in this region. The 
only significant backgrounds expected are from generic 
B and continuum events. The yields of the generic B 
decays are fixed to values expected from MC simulation 
while the continuum and non-resonant yields are allowed 
to float. There are 1100 data events in the selected Dalitz 
region and the fit yields —5.1 ± 7.6 non-resonant events. 
This is consistent with zero and the non-resonant contri- 
bution is therefore not considered as a background to our 
signal. 

An unbinned maximum likelihood fit to the variables 
Toes, AE, Ann and used to extract the total 

number of signal B^ — > p^ir^ and continuum background 
events and their respective charge asymmetries. The like- 
lihood for the selected sample is given by the product of 
the probability density functions (PDF) for each individ- 
ual candidate, multiplied by the Poisson factor: 



C=^^e-^' (iV')^n^- 



where A'' and N' are the number of observed and expected 
events, respectively. The PDF Vi for a given event z is a 
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sum of the signal and background terms: 



Sig 



where Qi is the charge of the pion in the event, 
^Sig(^Bkg>j ASig(^Bkg^ ^j^g yjgj^ asymmetry 

for signal (background) component j, respectively. The 
fractions of true signal (/sig), SCF signal (/scf), and 
wrong-charge mistag events (/mis) are fixed to the num- 
bers obtained from MC simulations. The j individual 
background terms comprise continuum, b c decays, 
p* , and seventeen other exclusive charmless B meson de- 
cay modes. Signal and continuum yields are allowed to 
float in the fit, with the generic B yields fixed to values 
expected from MC simulation. The PDF for each compo- 
nent, in turn, is the product of the PDFs for each of the fit 
input variables, V = V{mEs, ^E)V{At^N)V{m.,rTr)- E)ue 
to correlations between Ai? and mES: the 'P{mEs, ^E) 
for signal and all background from B meson decays are 
described by two-dimensional non-parametric PDFs [l^l 
obtained from MC events. For continuum background, 
V{mESy ^E) is the product of two one-dimensional non- 
parametric PDFs; TTiEs is well described by an empirical 
phase-space threshold function [T5| and AE is param- 
eterized with a second degree polynomial. The parame- 
ters of the continuum PDFs are allowed to float in the 
fit except for the endpoint of the empirical phase-space 
threshold function which is fixed at 5.29GeV/c^. Ann is 
described by the product of an exponential and a polyno- 
mial function for continuum background and by a Gaus- 
sian with a power-law tail on one side for all other 
modes. For 'P{mTr-n)^ one-dimensional non-parametric 
PDFs obtained from MC events are used to describe all 
modes except the signal mode itself, which is described by 
a non-relativistic Breit-Wigner line-shape. The param- 
eters for this PDF arc held fixed to the MC values and 
varied within errors to estimate systematic uncertainties. 

A number of cross checks confirm that the fit is unbi- 
ased. Using a double Gaussian PDF instead of a Breit- 
Wigner or omitting TOtttt altogether as a fit variable has 
no significant effect on the measured branching fraction. 
In 1000 MC pseudo-experiments, we use the maximum 
likelihood fit to extract the yields and asymmetries. The 
distributions for each component are generated from the 
component's PDF, giving values for the fit variables mES, 
A£^, Ann and TOtth- The expected number of events is 
calculated from the branching fraction and efficiency for 
each individual mode. The generated number of events 
for each fit component is determined by varying the ex- 
pected number according to a Poisson distribution. The 
test is repeated using samples with different asymmetry 
values. We repeat these MC studies using fully simu- 



lated signal B^ p^n^ events instead of generating the 
signal component from the PDFs. This verifies that the 
signal component is correctly modeled, including corre- 
lations between the fit variables. As another cross check 
we compare the distribution of the helicity angle ^hcI be- 
tween the momenta of the charged pion and the B meson 
in the p rest frame in data with that modeled in MC sam- 
ples for a variety of selection criteria. To investigate the 
possible effects of interference, we repeat the analysis ex- 
cluding events where both m^±^o combinations were in 
the range 0.55 to 0.95 GeV/c^; the branching fraction 
decreases by 0.1%. k 

TABLE I; Summary of the systematic uncertainties. 



Source 



Absolute uncertainties on yields 



^'"''^ (Events) 



B background normalization 


+ 6.9 
- 7.2 


PDF shapes 


+ 4.7 
- 4.2 


SCF fraction 


±12.2 


Mistag fraction 


± 2.0 


AE shift 


± 2.6 


Total 


±15 


Relative uncertainties on B{B^ - 




Source 


^Syst.(%) 


Efficiency estimation 


±7.3 


B counting 


±1.1 


Total 


±7.4 


Uncertainties on Acp 




Source 


'^Syst. 


Background normalization 


±0.006 


Background asymmetry 


±0.024 


PDF shapes 


±0.001 


Total 


±0.02 



Individual contributions to the systematic uncertainty 
are summarized in Table HI For each contributing exclu- 
sive B meson decay mode, we vary the number of events 
in the fit by its measured uncertainty, or by ±100% if 
derived from an upper limit. For the b —> c component, 
we fix the rate based on the number calculated from MC 
samples and vary the amount based on the statistical un- 
certainty on this number. The shifts in the fitted yields 
are calculated for each mode in turn and then added in 
quadrature to find the total systematic effect. To take 
into account the variation of the two-dimensional non- 
parametric PDFs used for AE and toes, we smear the 
MC-generated distributions from which the PDFs are de- 
rived. This is effectively done by varying the kernel band- 
width up to twice its original value. For m.n-ir and 
Ann, the parameterizations determined from fits to MC 
events are varied by one standard deviation. The system- 
atic uncertainties are determined using the altered PDFs 
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and fitting to the final data sample. The overall shifts in 
the central value are taken as the size of the systematic 
uncertainty. We vary the SCF fraction by a conservative 
estimate of its relative uncertainty (±10%) and assign 
the shift in the fitted number of signal events as the sys- 
tematic uncertainty of the SCF fraction. To account for 
differences in the neutral particle reconstruction between 
data and MC simulation, the signal PDF distribution in 
/S.E is offset by ±5 MeV and the data are then refit. The 
larger of the two shifts in the central value of the yield is 
2.6 events, which is taken as the systematic uncertainty 
for this effect. 




7l7t mass (GeV/c^) 



FIG. 1: Likelihood projection plots for the four fit variables, 
(a) mBS, (b) Ai?, (c) Ann, and (d) m^^. Each plot shows 
the total PDF (solid line), total background (dotted line), 
continuum contribution (dotted-dashed line), and the signal 
component (dashed line). 

Corrections to the 7r° energy resolution and efficiency, 
determined using various data control samples, add a sys- 
tematic uncertainty of 7.2%. A relative systematic un- 
certainty of 1% is assumed for the pion identification. A 
relative systematic uncertainty of 0.8% on the efficiency 
for a single charged track is applied. Adding all the above 
contributions in quadrature gives a relative systematic 
uncertainty on the branching fraction of 7.3%. Another 
contribution of 1.1% comes from the uncertainty on the 
total number of B events. 



To calculate the effects of systematic shifts in the 
charge asymmetries of background modes, the asymme- 
try of each mode is varied by its measured uncertainty. 
For contributions with no asymmetry measurement, we 
assume zero asymmetry and assign an uncertainty of 
20%, motivated by the largest charge asymmetry mea- 
sured in any mode so far |15| . The individual shifts are 
then added in quadrature to find the total systematic 
uncertainty. In addition, the effect of altering the nor- 
malizations of the B backgrounds affects the fitted asym- 
metry. The size of the shift on the fitted Acp is taken 
as the size of the systematic uncertainty. 

The central value of the signal yield from the max- 
imum likelihood fit is 365 ± 49 events, with a back- 
ground of 44840 ±217 continuum events and an ex- 
pected background of 842 ± 34 events from other B 
decays. Fig. [1] shows the distributions of rriESj 
Ann and m,rir- The plots are enhanced in signal by se- 
lecting only those events which exceed a threshold of 
0.1 (0.05 for Ann) for the likehhood ratio [16] R ^ 
(^Sig-pSig)/(^Sig-pSig _^ ^^^Bkg^Bkg^^ where N are 

the central values of the yields from the fit and V are 
the PDFs with the projected variable integrated out. 
This threshold is optimized by maximizing the ratio 
S = (TVSig eSig)/^7VSig gSig + ^.TvBkg ^Bkg ^j^^^.^ ^ 

the efficiencies after the threshold is applied. The PDF 
components are then scaled by the appropriate e. The 
efficiencies for the likelihood ratios vary for each vari- 
able and result in a different number of events in each 
projection. Compared against the null hypothesis, the 
statistical significance ^—2 \ii{CMuii I Cm ax) of the sig- 
nal yield amounts to 8.7 standard deviations. We ob- 
tain B{B^ p^n°) = [10.2 ± 1.4 ± 0.9] x 10-^ and 
Acp = —0.01 ± 0.13 ± 0.02, where the first error is sta- 
tistical and the second error systematic. The measure- 
ments are consistent with previous results [5] and provide 
improved constraints for the determination of the angle 
a from B —i- pi: decays. 
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